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Abstract. We use probability and several complex variable theory to study the density 
of complex zeros of a system of real random polynomials in several variables. We use the 
Poincare-Lelong formula to show that the density of complex zeros of a random polynomial 
system with real coefficients rapidly approaches the density of complex zeros in the complex 
coefficients case. We also show that the behavior the scaled density of complex zeros near 
of the system of real random polynomials is different in the m > 2 case than in the 
m = 1 case: the density goes to infinity instead of tending linearly to zero. 
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1. Introduction 



The density of real (resp. complex) zeros of random polynomials in one and several vari- 
ables with real (resp. complex) Gaussian coefficients has been studied by many. Kac [Kac48] 
and Rice [Ric54] independently found the density of zeros of a random polynomial with real 
standard Gaussian coefficients. Bogomolny, Bohigas, and Leboeuf ([BBL92], [BBL96]) and 
Hannay [Han96] have results on the density of (and correlations between) zeros of random 
polynomials with complex Gaussian coefficients. Edelman and Kostlan [EK95] generalize 
the results for density of real (resp. complex) zeros to systems of independent random func- 
tions in several variables when the coefficients are real (resp. complex) Gaussian random 
variables. 

In one variable, Shepp and Vanderbei [SV95], Ibragimov and Zeitouni [IZ97], and Prosen 
[Pro96] have studied complex zeros of real polynomials. Shepp and Vanderbei extended 
Kac's formula for the density of zeros of polynomials in one real variable, in the case where 
the coefficients are standard real Gaussian coefficients, to include both real and non-real 
zeros of those same polynomials. Ibragimov and Zeitouni studied the density of zeros of 
random polynomials with i.i.d coefficients (which are not necessarily Gaussian). Prosen 
followed Hannay's approach and found both an unsealed and a scaled density formula for 
the complex zeros of a random polynomial with independent real Gaussian coefficients. 

One consequence of Prosen's unsealed density formula is that, away from the real line, the 
density of complex zeros of an SO (2) random polynomial (which is the polynomial given by 



j=0 ^-^ ^ 

where aj is a real standard Gaussian random variable), rapidly approaches the density of 
complex zeros of a random SU(2) polynomial (which is the polynomial given by 



where Cj is a complex standard Gaussian random variable), as N, the degree of the polyno- 
mial, goes to infinity. Figure 1 illustrates this fact. 

Using the Poincare-Lelong formula, we show this convergence, recovering Prosen's single 
variable result [Pro96], and we show the convergence to be exponential. In Theorem 1 we 
generalize this result to the density of zeros of a random SO(m -|- 1) polynomial system in 
m variables (defined below). This generalized result is illustrated in Figure 2 for the special 
case m = 2. In [Mac09], the author uses a different method to generalize this result further 
to the density of critical points of a random SO(m + 1) polynomial in m variables. 

In one variable, Prosen also showed that, for every A^, the density of zeros tends linearly 
towards zero as we approach the real line. This result can be seen in Figure 1. In several 
variables, the behavior near \\y\\ =0 (i.e., near W^) is much different. See Figure 2. The 
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same can be said for the scaled density. In Theorem 2, we show that near M™, the behavior 
of the scaled density of complex zeros of the system in Theorem 1 is different in the m > 2 
case than the behavior of the scaled density in the m = 1 case that was shown by Prosen: 
for m > 2, the density goes to infinity instead of tending linearly to zero. Figure 3 illustrates 
the scaled density for the SO(m + 1) polynomial for m = 1 and m = 2. 

1.1. Density of zeros. Consider /im.Af = {fi,N, ■■■ , fm,N) '■ C'" — )■ C"^, where fq^N is an 
SO(m + 1) polynomial of the form 




where the Cj's are independent complex random variables, where the random vector {cj} e 
C^^,D^ = (^+'") has associated measure d'y, and where we are using standard multi-index 
notation. We show that for these m independent functions in m variables, the density of 
complex zeros in the real coefficients case rapidly approaches the density in the complex 
coefficients case as the degree of the polynomials gets large. In fact, we show that the 
convergence is exponential. Figure 2 illustrates this convergence for the case m = 2. 

Note that we have normalized these density functions so that the density in the complex 
coefficients case is the constant function 1, and that, because of invariance properties of 
the SO(m -|- 1) polynomial, it is enough to show the graph of the density function for z = 
{iy, 0) G C", < ?/ < 1. 

1.25- 




y~ 

Figure 1. The density of complex zeros of a random S0(2) polynomial for 
increasing values of A^. Because of symmetries, it is sufficient to plot the 
density along the imaginary axis for < y < 1. Here we have normalized so 
that the density of zeros of a random SU(2) polynomial is the constant function 
1. 
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More formally, let 




Figure 2. The density of complex zeros of the SO(m + 1) polynomial for 
increasing values of N , in the m = 2 case. Because of the invariance properties 
of the S0(m + 1) polynomial, it is sufficient to plot the density for z = {iy, 0) G 
with < y < 1. We have normalized so that the density of zeros of the 
random SU(m + 1) polynomial is the constant function 1. 



2.0—1 2.0 



1.5— 1.5 




Figure 3. The scaled density of zeros for m = 1 (left) and m = 2 (right). 
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where c G C'^'^, and S^dj^t is the delta measure on M-^^ c C'^^. Here d'jcx corresponds to the 
standard complex Gaussian coefficients case, where we are considering the random SU(m + l) 
polynomial 



|J|=0 ^ ^ 



where the Cj's are standard complex Gaussian random variables, and d'y^eai corresponds to 
the standard real Gaussian coefficients case, where we have the random S0(m+1) polynomial 

/ M\ 1/2 _N_ / j^j^ 1/2 ^ 



Z 



\J\=0 ^ ' \J\=0 

where Cj = Oj + iO is a standard real Gaussian random variable. Let E.y{-) denote the 
expectation with respect to 7; or, in other words, integration over C^^ with respect to d'j. 
Let Zfi^ = denote the distribution corresponding to the zeros of hm,N{z). Then 

E^{Zh^ ^(2)) denotes the density of the zeros of h with respect to the measure d'j. We now 
formally state the ffist result: 

Theorem 1. 

for all z G C™\]R™, where is a positive constant that depends continuously on z. The 
explicit formula for is 

A, = - lo£ 



1 + z- z 



1 + Ikll 

Also, for compact sets K C C™\]R™, the density converges uniformly with an error term of 
0{e~^^^), where Xk is a constant that depends only on K . 

Note that for z G C"*\]R™', the argument of the log is less than 1, and is positive. The 
formula for E^^^{Zhj^(^z)) is a special case of a result in [EK95], and is a very simple function: 

F (7 1 - 1 

The formula for E^^^^^{Zh^^z)) is very complicated, but, by Theorem 1, we know that 
E^^^^^{ZhJ^{z)) equals a very simple function, E^^^{Zh^(^z)), plus some exponentially small 
term. 

We prove this result using the Poincare-Lelong formula, which is similar to that which 
was used in [BSZOOa], but has the added complication that the coefficients are real. The 
proof uses 2-point Szego kernel asymptotics, which still applies to the polynomials with real 
coefficients because we are viewing them as functions of complex variables. We also use the 
fact that the /q^7v(^)'s are independent, which is a major difference from the critical points 
case considered in [Mac09]. 

Shiffman and Zelditch [SZ99] and Bleher, Shiffman, and Zelditch ([BSZOOa], and [BSZOOb]) 
have generalized many results about random polynomials on C"^ and to complex mani- 
folds, and they have several results related to the statistics of zeros of a random holomorphic 
section of a power of a line bundle over a complex manifold. In particular, in [BSZOOa], 
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the authors use the Poincare-Lelong formula to find a formula for the density of zeros and 
correlations between zeros. Edelman and Kostlan used a similar approach in [EK95] to get 
a result like E^^^{Zhj^(^z)) but for systems of more general complex functions. 

In [DSZ04], Douglas, Shiffman, and Zelditch look at the critical points of a holomorphic 
section of a line bundle over a complex manifold, motivated by applications in string theory. 
They use a generalized Kac-Rice formula to find statistics of these critical points, namely 
the density of critical points and correlations between critical points. In [Mac09], we study 
complex critical points of a random polynomial with real coefficients and generalize the result 
in Theorem 1 to the density of critical points of a SO(m + 1) polynomial. 

1.2. Behavior of the scaled zero density near M™. Consider hm,N{z) as above. Note 
that the behavior of the density function in the m = 1 case (Figure 1) and the m = 2 case 
(Figure 2) differ greatly. Consider the scaling limit of the density. 




which will help us understand the behavior of the density function in a region around 
that is shrinking at a rate of We can show that K^^^^^[z) depends only on ?/ = Imz, so 
we can write the scaled density as K^^^^{y). 

Figure 3 illustrates the behavior of (y) near for m = 1 and m = 2. Note that 
for m = 2, -f^^^^, (y) does not tend linearly towards zero as in the m = 1 case, but instead it 
tends to infinity. We prove the following: 

Theorem 2. For y near 0, 



We will see that difference in the m = 1 and m > 2 cases boils down to the fact that 



Finally, after working mostly on C\M and C™'\]R™', we give a weak limit for E.y^^^^{Zfj^{z)) 
in the m = 1 case. We show that j^E^^^^^{Zfj^{z)) = j^E^^^{Zfj^{z)) + 0{N~^), weakly on 
K C C We stress that K could contain some points in M, whereas our strong convergence 
result excludes points in M. 



KZJy) = oi\y\) 



m = 1, 





for y e M\{0}, 




2. Proof of Theorem 1 for m 



1 



Consider the real random polynomial 



N 
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where the a^s are real independent Gaussian random variable with mean and variance 
(^) . Alternatively, one often writes 



Z 



£=0 

where ai is a standard real Gaussian random variable. Instead, we choose to think of the 
random polynomial 

^\ r. / 

e=o 

where q is a more general complex random variable with associated measure d'j. We then 
consider two special cases 



d^real = 5s^e-\'^\'dc,CeC''+\ 

where Ss is the delta function on 5* C C^"^^, the set of points c = a + z6 G C^^^ where 
6 = e M^"*"^. Here djcx corresponds to the standard complex Gaussian coefficients case, 
where we are considering 

z' 

- V /' / 

e=o 

where the q's are standard complex Gaussian random variables, and d'jreai corresponds to 
the standard real Gaussian coefficients case, where we have 



fNiz) = 2^cA z =2^a, z , 

e=o ^ e=o ^ ^ 

where q = + iO is a standard real Gaussian random variable. We let E^{-) denote 
expectation with respect to dj. 

The goal of this section is to prove Theorem 1 in the m = 1 case: 

Proposition 2.1 (Theorem 1 for m = 1). We can write 

E,..AZfA^)) = E.^^iZjJz)) + Er^iz), 

where 

N 1 



E,..iZf(z)) 



I2^2 



TT (1 - , _ 

and where Epf{z) = 0(e~'^^^), for all z G C\]R. Eere \z = — log | x+jfp| ^■^ ^ positive constant 
that depends continuously on z. 

Using the Poincare-Lelong formula, we write E^^^^^{Zf^{z)) = E^^^{Zf^{z)) + E]sr{z), and 
we then aim to find an explicit formula for this error term. Writing fN{z) = a ■ F]sf{z) and 
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un{z) 



, we can write 



ll^ivWII 



En{z) dx a dy = E. 




71 



I 



dd\og \a ■ um{z)\ 



) 



Note that E^^^^^{-) denotes an integral over all of with respect to d'jreai and is fairly 
complicated integral. However, also note that d'jreai, the real standard Gaussian measure, is 
rotationally invariant. We use this fact and perform 2 real orthogonal changes of variables 
in order to simplify this integral over and write it as an integral over M^: 

Lemma 2.2 (Performing real rotations). 



for some functions r = r^^z), s = sn{z), and t = t^iz) for which we give formulas within 
the proof of the lemma. 

We make another change of variables, this time switching (oq, ai) to polar coordinates 
{p,0). We can easily integrate with respect to p, further simplifying our error term to an 
integral with respect to 9 over [0, 27i]. We apply Jensen's formula to what is left to evaluate 
the integral and get an explicit formula for En{z): 

Lemma 2.3 (Evaluating the integral), //r^ + + = 1 and r,s,t > 0, we have 



Lemma 2.4 (An exact formula for, and asymptotics for, the error term). We show that, by 
Lemma 2.2 and Lemma 2.3, and after some simplification of2rt, we have 



It follows that En{z) = 0(e ^''^),z E C\]R, where = — log | ^^^Ip \ is a positive constant 
depending continuously on z. 

This result gives us Theorem 1 for the one variable case. 

In subsections 1 through 5, we prove the aforementioned results. The approach we use 
is similar to that described in [BSZOOa], where they find the limit of the pair correlations 
of zeros of random holomorphic sections of powers of a line bundle of a complex manifold. 
While we only deal with density of zeros in this section, the condition that the coefficients 
aj are real causes the method in [BSZOOa] to be useful. 

2.1. Proof of Proposition 2.1 - Theorem 1 for m = 1. We write a = (ag, ...,a7v) and 
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SO that /at = a-F^. By the Poincare-Lelong formula, the density of the zeros of /, E^^^^^{Zf^), 
satisfies 

E,_,iZf,) dxAdy = E,_,C-dd\og = E,_,i-dd\og \a ■ F^\). (1) 

TT 71 

We write F^i^z) = \ \F]\f{z)\\u]\f{z), where U]\f{z) is a unit vector, and (1) becomes 
E,_,i-dd\og\\F,,iz)\\) + E,_,i-dd\og\a.u^iz)\) 

71 TT 

From [BSZOOa] we can see that 

E^..iZf^iz)) dx Ady = -dd\og\\Fr,{z)\\ = E^_^i-ddlog\\Fr,{z)\\) dxAdy, 

71 IT 

SO we have 

We now prove the Lemmas, which deal with the error term Em{z)^ and return to the proof 
of Proposition 2.1 in Section 2.5. 

2.2. Proof of Lemma 2.2 - Real Rotations. As shown in [BSZOOa], when aj is a standard 
complex Gaussian random variable, this error term Ejq{z) is zero for all N (not just as 
— )■ oo). Because of the SU(2)-invariance of the standard complex Gaussian measure, one 
can perform a unitary change of variables so that u becomes (1,0, ...,0) and the integral 
j^N+i ddlog \a ■ u\d'^cx{.ci) becomes a single integral that evaluates to 0: 

/ (991og|a- (l,0,...,0)|d7c^(a) = / 991og |ao|rf7cx(ao) = 0. 

In the case where aj is real, the second term is not zero for all A^. Because only real 
rotations can be performed, u can not be rotated to (1,0, ...,0), giving a single integral. 
But we can still use the rotational invariance of real Gaussian measures to obtain a double 
integral over M^, which is a little more manageable than the integral over M^"*"^. 

Let u = Reu + ilmu = ( Re Ui, Reu^) + i( Imwi, Iirun). Note that u, Re u, and 
Imu depend on z and A^ but we frequently omit these arguments for convenience. Since 
we need to do real rotations, the real and imaginary parts of u must be rotated the same. 
Therefore, as mentioned, we can not rotate u to (1, 0, 0). However, we can rotate 
so that either the real part or the imaginary part of u is of the form (r, 0,...,0), where 
r = r^lz) is some (non-zero) constant less than 1. So we choose to perform a (real) rotation 
of oq, ai, ajv so that 

u = Reu + ilmu = (r, 0, 0) + i{lmui, ImuN). 

Then one can perform a rotation of the Oi, ajv variables so that Reu is unaffected and u 
becomes 

{rN{z), 0, 0) + i{sN{z),tN{z), 0, 0) 
= {rNiz) + isN{z),itN{z), 0, 0). 
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Note that since m is a unit vector, and rotations preserve length, r, s, and t have the condition 
+ + = 1. Note also that r, s, and t all depend on z and but we frequently omit 
these. We are now concerned with the limit of the simpler integral. 



TT 

i 



En{z) =- I dd\og\{ao,ai,...,aN)-ir + is,it,0,...,0)\d-freai{a) 
(99 log \ao{r + is) + ai{it) \ dii{ao)dfi{ai) 
= -dd I \og\ao{r + is) + ai{it)\—e''^''o+^^^^/^daodai. 

2.2.1. Formula for r. . First, we know that since u{z) = and since the length of Rew 

doesn't change from a rotation, we can write 

r^ = ||Re«|r = ||Ret.|r = L_^ = _ + _Re (^^^j = 3 + Ole"^^"'), 

where e"'^'' = [ i^\f\2 \ < I, z E C\]R, Note that we are assuming Imz 7^ 0, and note that 
Xz = - log 



l+z- 



> 1. 



2.2.2. Formula for s.. Next, we have the relationship Rew ■ Imw = rs, so since the angle 
between Rew and Im-u doesn't change under a rotation, we have 

Re u ■ Im u Re m ■ Im u \^\ ' jj^ Re F ■ Im F 1 



IFIP r 



u ri + z'V 1 



2lm(^^^) .- = 0{e'''^),zeC\R. 

2.2.3. Formula for t.. Since + + = 1, we have t easily: 

[tN{z)]' = 1 - [r^(^)]2 - [s;v W]' = 1 + 0(e-^^^), ^ G C\M. 

Also, we note that r, s, and t converge uniformly on compact sets K C C\]R, where the 
error term is 0(e~^^^), for some universal constant A (independent of and z). Finally, we 
note that all derivatives of r,s, and t are 0(e~'*'^^), and all derivatives converge uniformly 
on compact sets K C C\]R as well. 



2.3. Proof of Lemma 2.3 - Evaluating the integral. We now prepare to use Jensen's 
formula to evaluate the integral 

-dd I log|ao(r + is) + ai(2t)|— e-(''"+''?)/2^aoC^ai- 
71 Jjj2 2tt 



10 



BRIAN MACDONALD 



2.3.1. Switch to polar coordinates. First, we switch to polar coordinates, so the integral 
becomes 

'"27r POO 



27r2 



dd / log \{p sin 9) {r + is) + {p cose) {it)\e-^ pdpde. 

J 0=0 Jp=0 

We may factor out a p from the argument of the log and get 

'•2tt poo 



dd 



(log p + log I (sin 9) (r + is) + (cos 9) (it) \ ) e"^ pdpd9 



Jp=0 



Since 



2tt poo 



logpe pdpd9 



J p=Q 

doesn't depend on it gets killed by 55, so we are left with 

p2tt poo 

dd j / log|(sin^)(r + + (cos^)(it)|e-'''/Vc^pd^. 

J 6=0 Jp=0 

The log term doesn't depend on p, so we may pull that term outside the integral, and 
integrate with respect to p to get 



27r2 



27r2 
i 

i 

i 



p2tt 

dd log I (sin 61) (r + is) + (cos 61) (it) I 

Je=o iJp=o 

/■27r 

dd / log I (sin 9) (r + is) + (cos 9) (it) \ 
Je=o 

/■27r 

dd / log I (sin 9) (r + is) + (cos 9) {U) \ [1] d9 
Je=o 

p2it 

dd / log I (sin 9) (r + is) + (cos 9) {U) \ d9. 
'e=o 



-pV2 



pdp 



d9 



pV2 



d9 



2.3.2. Jensen's Formula. Using the fact that cos ^ = |(e*^ + e '^) and sin ^ = |t (e*^ — e *^)., 



we can write 



27r2 



27r 



1 



dd / \og-\{e''' -e-'''){-i){r + is) + {e'' + e-'yt\d9. 



e=o 



We can bring out a log |, and since 95 log | = 0, we have 



27r2 



27r 



55 / log \(e'' - e-*')(-ir + s) + {e'" + e'^ytl d9. 



e=o 



We can factor out e , and since log |e | = 0, we get 

f27r 



27r2 
i 

2^ 



55 



55 



=0 



log 1(6*2" - l){-ir + s) + (e^2" + l)it\ d9 



2n 



log |(s + i{t - r))e'^^ + i-s + i{t + r))\ d9. 



9=0 
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We can now use Jensen's formula to evaluate the inner integral. Recall that Jensen's formula 
states that, assuming 0(0) 7^ 0, and is non-zero on dD{0, 1), then 

'^og|0(e^^)|rf^ = log 10(0)1+ V log^ 



1 

2^ 



In our case, 



so that we have 



[w 



) = {s + i{t- r))w^ + (-S + i{t + r)), 



0(0) 



{s + i{t - r))e^^^ + {-s + i{t + r)) 
—s + i{t + r), and 

—s + i{t + r) 







s + i{t-r) 

Note that since |0(O) p = | - s + z(t + r) p = + (t + r)^ = + + 2rt + = 1 + 2rt, and 
r and t are non-negative (by construction), 0(0) 7^ 0. 

We show that |wj(2;)p > 1, for all z, implying that |wj(2;)| > 1 and that all the zeros Wj{z) 
of are outside the unit disk for every z. We have 

-s + i{t + r) 2 s"^ + + 2rt + 



\Wj{z)\' 



s + i{t — r) 



+ t'^ - 2rt + 



> 1 



for z E C since r and t are non-negative by construction. Note that the only time rt is zero 
is when 2; G M. In this case, r^iz) = 1 and Sn{z) = tN^z) = for all A^, and |wj(2;)| = 1. 
So since all of the zeros of are outside the unit disk, we have for 2 G C, 

r log \<l){e'')\de = log 10(0)1 =\og\-s + t{t + r)\ 

= hog\- s + z{t + r)\^ = ^log(l + 2rt). 



or 



1 

2^ 



log \ao{r + is) + ai{it)\e "'(iaoc/ai = ^ log(l + 2rt), 



the result of the lemma. 

2.4. Proof of Lemma 2.4 - An exact formula for the error term. From the proof of 
Lemma 2.3, we have 



E]\f{z) dx A dy 



27r2 
i 



dd I log|0(e^^)|de 



--^9a(27r)^log(l + 2rt) 



1 92 

—dd\ogil + 2rt) = - — 
2tt tt ozoz 



- log(l + 2rt) dx Ady,z e C\M. 



After some simplification of 2rt we have 

1 92 , 



En z] 



7T dzdz 







(1- 












\z\^)^ 


1 






(1 + 





, z G C\M. 
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2.5. Finishing Proof of Proposition 2.1 - Theorem 1 for m = 1. Since we found that 

r = rj^iz) = ^ + 0{e-^^^),t = t^iz) = 0(6"^-^), and s = Sn{z) = 0(e-^^^) and 

all derivatives (in particular, the first and second derivatives) of r, s, and t are 0(e~'^"'^), we 
can say that, from the proof of Lemma 2.4, 

Ej^{z) = 0{e~^'^),ze C\R. 

The proof follows directly from this result, combined with Lemma ?? and Lemma 2.4. 

3. Proof of Theorem 1 for m > 2 

In this section we are concerned with the zeros of hm^N = {fi,N, fm,N) '■ C'" — )■ C™, 
where fg^^ is a polynomial of the form 




where Oj is a real standard Gaussian random variable, and where we use the following 
multi-index notation: 

Z (^Zi , . . . , Zffi^ 

\J\=ji^ h 

a''j = al ■ eR 
/N\ _ f N \ _ N\ 

\j J {N - ji- ... - j^y.jii ... jj. 

^ ^1 •••^m • 



Instead, we choose to think of the random polynomials 




where the Cj's are independent complex random variables, where the random vector {cj} G 
Dn,Dn = (^^'") has associated measure d'y for each q. We then consider two special cases 

7t" \ m J 

where 5s is the delta function on S* C C^^, the set of points c = a + ib G C^'^ where 
b = E M^^. Here djcx corresponds to the standard complex Gaussian coefficients case, 
where we are considering 
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where the Cj's are standard complex Gaussian random variables, and d'jreai corresponds to 
the standard real Gaussian coefficients case, where we have 

/^^^ 1/2 /^^^ 1/2 ^ 



z 



\J\=0 ^ ' |J|=0 

where Cj = Oj + is a standard real Gaussian random variable. We let denote 
expectation with respect to d'^reai and E^^^{-) denote expectation with respect to djcx- 
The goal of this chapter is to prove Theorem 1 in the m > 1 case. 

Theorem 1. 

for all z G C™\]R™, where is a positive constant that depends continuously on z. The 
explicit formula for is 

Xz = - lo£ 



1 + z- z 



12 



Also, for compact sets K C C™\]R™', the density converges uniformly with an error term of 
0(e~'^^^), where Xk is a constant that depends only on K . 

So at any point away from M"*, the expected density of zeros in the real coefficients case 
approaches the expected density of zeros in the complex coefficients case as N gets large. 
Using the Poincare-Lelong formula, we can write 

where 



and where En{z) is some "error term." We then find an explicit formula for this error term. 
Writing fq,N{z) = a ■ Fq^N{z) and Uq^N{z) = |j^'^^fy||; "we can write En{z) dx A dy as a sum 
of 2™ — 1 terms, each of which contains a factor of the form 

^w. (zi^n^:^log|«-«^WI) (2) 



Note that E^^^^^{-) denotes an integral over all of M^^, with respect to d'jreai, and is a fairly 
complicated integral. However, like in the one variable case, we note that d'jreai is rotationally 
invariant, and we use this fact and perform 2 real orthogonal changes of variables in order 
to simplify this integral over and write it as an integral over M^: 

Lemma 3.2 (Lemma 2.2 for m > 1). The factor (2) equals 

1 C 1 

- / log |ao(r + zs) + ai(zt)| — e~^'^o+"-^^)/'^ daodai 

TT J^2 271 

for some functions r = r^lz), s = Sn{z), and t = tN{z) for which we give formulas within 
the proof of the lemma. 

Then we easily get 
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Lemma 3.3 (Lemma 2.3 for m > 1). . //r^ + + = 1 and r,s,t > 0, we have that 

If 1 

— / log|ao(r + is) + ai(i)f:)|e~^"o+"i^/^rfaoC?ai = -log(l + 2rt). 
2n J^2 2 

This means that the factors (2) are of the form 

1 

TT OZjOZk 

each of which is 0{e~^^^) for some positive constant Xz depending continuously on z. Using 
an argument that rehes on the independence of the fq,N{,z)'s, we get the following: 

Lemma 3.4 (Lemma 2.4 for m > 1). We give an exact formula for Ei^{z), and show that 
it goes to zero rapidly, i.e., 

En{z) = 0(e-^^^), z E C^XW^, 

where Xz is a positive constant depending continuously on z. We give an explicit formula for 
Xz in the proof of the theorem. 

This lemma gives us our first main theorem, as shown in Section 3.5. 
3.1. Proof of Theorem 1. We begin by writing 
«' = M „,..., <i'j,...,a2 „,„)eK°" 



where = m™)' ^^^^ write fq^j\[ = o-'^ ■ Fq,N = ci'^ ■ -^tv- 

By the Poincare-Lelong formula, we have 

EiZj.iz) A ... A Zfjz)) = Ei^ddhg A ... A ^ddlog \f^\') 



E 



2^ 



log ■ F|2 A ... A 99 log la*" ■ Fj^) 



which we can write more succinctly as 

m m 



-j E(^l\dd\og\a'' ■ F\ 

' q=l 



Writing F = Tr^||F|| and u = yr^, we can write (3) as 



2^ 

i 

2^ 



m m 



E(/\dd\og 



q=l 

m m 



e \ /\ {dd\og\\F\\^ + dd\og\a'^ ■ u\'^) 



q=l 



(3) 



(4) 
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Since 991og||F|p + ddlog\a'^ ■ is independent of for / 7^ q, then we can write the 
expected value of the wedge product as the wedge product of the expected values, and (4) 
becomes 



i 

2^ 



/\ E^dB log WFW^ + dd log \a'' ■ 

ddlog ■ 



9=1 

m m 



/\ (ddlog WFW^ + E 

q=l 



At this point, we could find the large limit; we have essentially reduced the m- variables 

case to showing that E ddlog \a'^ ■ is 0(e~^^^). If this term is indeed 0{e~''^^^), then all 

but one term in the wedge product goes to zero exponentially fast. Since we want an exact 
formula for the density of zeros, we delay the proof of the large limit, and we first work 
out the details of writing an exact formula more explicitly. From that formula, the large A^ 
limit and the scaling limit will follow easily. 
We write 



where 



Ei^N{z)du 
E2,N{z)dui :-- 



I 

2^ 



m m 



I 

i 

2^ 



Eyl\ ((9(91og||F||2 + 9(91og|a«-Mp) 

[Ei^n{z) + E2,n{z) + ... + E2rn^N{z)]duj 



E(ddlog\\F\\^ AddlogWFW'^ A ... ^ddlog\\F\\^ 
E(ddlog\a^ ■ u\ Addlog\\F\\'^ A ... Addlog\\F\\^ 



E2^,N{z)duj := ( — ] E(ddlog\a^ ■ u\ A ddlog\a^ ■ u\ A ... A ddlogla"^ ■ u\ 



We look at these 2™ terms and we claim that only the first term is non-zero in the limit. 
The first term is known: 

Ei^N{z)duj = E^^^ (Zh^) = ( ^ j E(^dd log I |F| |2 A dd log \\F\\^ A ... A ddlog \ \F\f^ 



1 



TT™ (1 + |zi|2 + ... + |z„ 



\2\m+l 



duo 



1 



duo 



We let the error term En{z) be the sum of the remaining terms: 

En{z) = E2^n{z) + ... + E2m^N{z). 

We now prove the lemmas, and continue the proof of Theorem 1 in Section 3.5. 
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3.2. Proof of Lemma 3.2. Consider the integral 

log ■ UNiz)\ dfi{a^). 

Note that the terms E2^n{z), -E2™,Af(^) each contain a factor of this form. Also note that 
this integral is in the same form as the error term in the one variable case (see Lemma 2.2), 
so we proceed in a similar manner. We rotate (ao-.-o, ao--oAf) and then (aio-.-o, ao--o7v) so 
that the integral becomes 

log \a^ ■ (r + is, it, 0, 0)| dfi{a^) 

log |ao...o('^ + is) + aio-oitl dfi{ao...o) rf/i(aio...o) 



where r = ri^{z),s = SAr(z), and t = t]^{z). We have similar formulae and large limits for 
r, s, and t: 



l + zl + --- + zl 

1 + |Zi|2 H \- \Zrr 



N 



+ 0{e-^^^),zE C 



Sn{z 



Im 



l + zl + --- + z 



1 + blP + --- + \Zr, 



N 



rNiz) 



[tr,iz)Y = 1 - [rNiz)Y - [sNiz)Y = - + 0(6-'^''), z e C"^\M™. 

Also, all derivatives of r, s and t are 0(e-^-^), z e C'"\M"'. Finally, we note that, r, s,t and 
their derivatives converge uniformly on compact sets K C C''"\]R'", where the error term is 
(9(^g-Aif Af-j ^Qj, some universal constant A (independent of A^ and z). 

3.3. Proof of Lemma 3.3. By Lemma 2.3, we have 

If 1 

— / log |ao-o('^ + is) + aio...ozt\ rf/i(ao...o) rf/i(aio...o) = 7, log(l + 2rt). 
Itt Jjj2 2 

3.4. Proof of Lemma 3.4. Consider, Eg the g-th term in En{z). This term is of the form 

Eq^N{z)duj = E{dd^\ A ... A 



where (pl j^{z) is either log ||F/v(2;)|| or log \a^ ■ U]^{z) \ for each £. For example, for E2^n{z) we 
have (pl = log \a^ ■ u\ and (pj = log ||F/v(-2)|| for 1 < ^ < A;. Writing out the wedge product 
we get 



En,Niz) 



E 



dz„[X)OZriX) 



92 



dz^{Tn)dzT-(^rn) ™ 
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where the sum is over all permutations a and r of {1, 2, m}, and where (—1)'^ denotes the 
sign associated to the permutation a. Since the sum is finite, we can write 



-0? 



or 



dZa{l)dZr{l) 



where E^^iz) := E 



n 



To simplify notation even more, let 



~ dz.<,,dz 



dza(e)dzr{e) 



so that we have 



EZAz) = E 



.e=i 



Now, note that (f)l{z) does not depend on all of a^,...,a^, but only depends at most 
on a^. (If = log ||FAr(z)||, then it doesn't depend on cither.) So because (p'^iz) is 

independent of for all i' ^ i, we will write this integral over M.^'^ x • • • x M'^^ as a product 
of integrals over M.^^ . To do this, we first let 

Lq = {i : (f)^ is of the form log ||F/v(2;)||} C {1, . . . 

Note that is of the form log ■un{z)\}, for all I ^ Lq. We can now we split the product 

to get 



E:m = E 



^E 



n 



n ^^'^K^) 



llDZlog\a'-UN{z)\ 



By the definition of expected value, we have 







/ 


[( 







n^riogii^iv( 



dii{a^)---dii{a"'), 



Note that the first product is independent of for all £ ^ Lq, and the second product is 
independent of all £ G L^, so we can write E^i^z) as 



llD'',^^log\\Fr,{z)\\dii{c 



Y\_ D'^'^ log \a^ ■ UN {z)\ diJi{a^) 
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The first product is also independent of for all £ G Lg, and since /j^d^v dii[a^) = 1, we have 
for the first integral 

Even more, the £-th factor in the second product depends only on i and is therefore inde- 
pendent of all i' ^ L not equal to i. So the integral of this product becomes a product of the 
integrals: 

Y\ / D"^'^ log \a^ ■ UNiz)\ dn(a^) 
and we can switch the derivatives and the integral to get 



Putting everything together we have 



l[DZ log \\Fm{z) 



Y\ [ \og\a' ■ un{z)\ dfi{a') 



Lemma 3.3 gives us 



l[D;nog\\Fr,{z) 



3.4.1. Exact formula. Further simplification gives 



n ^rl 



log(l + \ \z\ 



\2\N 



.2 IUII2 



2\N 



+ 

where we use the notation z ■ z = zf + ■ ■ ■ + z"^, = + ■ ■ ■ + |-2mP- Finally, 
and 

2™ 

En{z) = Y^EqA^)- 
q=2 

3.4.2. Limit formula. All derivatives of log||F|| are bounded. Next, r = r]\i(z) = + 

0{e-^-^), t = t^iz) = ^+ 0{e-^'^), and s = SNiz) = 0{e-^-^) and all derivatives (in 

particular, the first and second derivatives) of r, s, and t are 0(e~'^"=^) on C^yM™. So we can 
say that all second derivatives of log(l + 2rt) are 0{e~^^^) on C'"\M™'. This means that 
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Since this is true for each i, a, and r, we have 

2'" 

q=2 cr,T 

3.5. Finishing Proof of Theorem 1. Using our work in Section 3.1 with Lemma 3.4, we 
get our main result: 

4. Proof of Theorem 2 

We then wish to study the behavior of E^^^^^{Zf^{z)) near the real line. We define the 
scaling limit of E^{Z fj^{z)) to be 

K^{z) = lim ^EJZff^)) 



The scaling limit helps us understand the behavior of the density function in a region around 

1 



that is shrinking at a rate of -4^. We note that 



' IT 

and we find the scaling limit of the error term when m = 1: 



Lemma 4.1 (Scaling limit for the error term, m = 1] 



E^{z) = lim ^Em{^) = log I 1 + Wl - 



r2 



2 



, z G C\M. 



By setting z = x + iy, we can write 

lim ^E^i^) = ;M^log (l + Vl - e-^y') ,y^0. 

Then, by Proposition 2.1 and Lemma 4.1, we recover Prosen's scaled density equation, 
and our Theorem 2 for the special case m = 1: 

Proposition 4.2 (Equation (26) in [Pro96], and Theorem 2 for m = 1). We have 

K- 1 wv + i)e:!!! 

lreal\ > ^ ^ _ g-4j/2 -J 3/2 

Since K!^^^^{z) depends only on y, we can write K^^^^^{y), and we have the asymptotics 

KZJy) = o{\y\) 

for y near zero. 
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Since y = corresponds to the real line, this result tells us that the scaled density tends 
linearly toward as we approach the real line. This formula for K^^^^ {z) was given by Prosen 
as mentioned, but we find it here using Poincare-Lelong method that will be generalized for 
the m > 2 case. We give a formula for the scaling limit of the "error term" E]\f{z), when 
z G C™\]R™, which we denote Eoo{z), z G C™\]R™, and the scaling limit of the density, which 
we denote, K!^^^^{y). (Once again, the scaled density only depends on y = Im^;.) We show 

for m > 2 is different than the behavior for m = 1: 



that the behavior of K°° (y) near 



Theorem 2. For y near 0, 



yreal 



iy) = o{\y\), 
1 



\y\ 



m = 1 
m> 2. 



We also show that as | Imz| — )• oo, Eoo{z) — )■ 0, so that K^^^^{y) — > 0, as ?/ — )■ oo. In other 
words, the scaled density of zeros in the real coefficients case approaches the scaled density 
of zeros in the complex coefficients case as you move far away from M"^. 

Finally, after working mostly on C\]R and C™\]R™, we give a weak limit of the error term 
Ej^{z) on compact sets C C (which may include points in M): 

Proposition 4.3 (A weak limit, m = 1). 



1 ^ 



N 



En{z) dx a dy = 0{N ^), weakly on compact sets K C C, 



by which we mean that for any G C°°{K), 



N 

This means that 
1 



{En{z) dx A dy, (p{z)) 



1 

N 



En{z) (piz) dxAdy = 0{N~^). 



K 



-E 



Ireal ' 



N 

1 

7r(l + \z\^y 



Zf,{z)) + 0{N~') 

+ 0{N-^), weakly on C C. 



Note that K could contain some points in 
points in M. 



whereas the strong convergence result excludes 



4.1. Proof of Lemma 4.1 - Scaling limit of the error term. By the chain rule we have 
for any differentiable function f{z) 



dzdz 

So we have from Lemma 2.3 

1 92 



N 



dzdz 



/( 



N 



— Em 

N ^ 



N 



Nti dzdz 
1 d^ 



-log [l + 2rN{z)tN{z)] 



TT dzdz 



-log 



1 + 2r 



z . , z 



N 



N 



z G C\M, 
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and after some simplification we get 

1 92 
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-F ( — 



We now take the limit and get 



3 log 



1 + 



v 



1 



d + \^\2)N 



,z e C\M. 



1 92 



zlog 1 + Wl 



, z e C\M. 



Setting z = X + iy, we can write 

lim ^^^(^) = -^1^ log (1 + Vl - e-^y') ,y^0, 
and after simplification and adding to K^^^{z) = ^, we recover Prosen's result in [Pro96]: 



ireal ^ ' 



1 1 - (%2 + i)e- 
~ (l_e-V)3/2 



Ay'' 



4.2. Proof of Proposition 4.2 - Theorem 2 for m = 1. Using Lemma 4.1 we have the 
asymptotics 



^ZM = o{\y\) 



for y near zero. 



4.3. Proof of Theorem 2. From the proof of Lemma 3.4 we have 



n 



n^r'2'og|i + t/i 



(5) 
(6) 



where we use the notation z-z = z1^ h-z^, 

then z ■ z = \ \x\ 



2 = \zl\^^ h |-2mP- If we write z = x + iy, 

P + 2i(a; ■ ?/) + I |y| p. Since je^**^^'^^! = 1, we can write the second product as 

n^r^log (l + Vl-e-4NP) 



Since the first product is bounded (it is 1 if a{i) = T{i), for all i & L, and zero otherwise), 
and the second product goes to zero exponentially fast as \ \y\ \ 00, we have 

Eq;L{z) 0, as \ \y\\^ 00. 

Since this is true for each i, a, and r we have 

2™ 

5=2 a,T 
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and since 

(z)= lim (Z.^^))= lim ^ ( 7v(^) + ^iv(^) ) , 

we get 

TTl 

(z) — , as\\y\\^ oo 

We now look at the behavior of the second term in 5 near W^. We have the following 
asymptotics: 

1 iog(i + Vi - exp-ibiP) = Hall + Odlall^). 
^l,og(i + yr^^ip^) . + 0(1) ^ o(^). 

Since in the worst case the second term has m products, we have at worst 



n«?1iog(i+/r3^)=o(^ 



givmg us 



\y\ 



q=2 
= 0(l 



q=2 a,T 
1 



4.4. Proof of Proposition 4.3 - Weak limit. Let C C be a compact set. Note that 
unlike before, we are including points on the real line. We now show that jjE]^[z) dx A dy 
goes to weakly on K. More specifically, we show that for any G C'^{K), 

^{En{z) dx A dy, (t){z)) = j^J^ Em{z) (t){z) dx A dy = O(iV-i). 
Recall that E]\f{z) dx A dy = E{—ddlog \a ■ un{z)\). By the definition of the expectation 

TT 

of a distribution, we have 

{En{z) dx A dy, (f){z)) = ( E{-dd\og \a ■ un{z)\) , (f){z) 

V ^ 

= E ( —99 log \a ■ un{z) \ , (f){i 
By the definition of the derivative of a distribution, we have 

E ( —ddlog \a ■ un{z)\ , (piz) j = E (log\a ■ un{z)\ , —dd(f){z] 
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and by the definition of a distribution, we can write tliis as 

E (^j log \a ■ UN {z)\ ^dd4){z] 

Recall tliat E denotes expectation witli respect to tlie Gaussian measure d^reai- We tlien 
liave by definition of expected value tliat tliis equals 

log |a ■ un{z) \ -dd(j){z) ] d'Jrealia)- 

K TT y 

Since the integrand is bounded, and since (^(z) does not depend on a, we can switch the 
order of the integrals and get 

/ log |a ■ UN{z)\d'yreai{a) ) -dd(j){z). 

K \JW* J ^ 

Recall that by Lemma 2.3 above we have that the inner integral is |log(l + 2rt), so we 
have 

log \a ■ UN{z)\d'yreai{.a) ) -dd(f){z) = / 7^ log(l + 2rt) dd(f){z) 

Recall also that rN{z) and tN{z) are both non- negative by construction, and both are 
bounded by 1 since + + = 1. Both of these conditions are true even on the real line, 
where rjv = and In = ^ for all A^. This implies the crude estimate 1 < (1 + 2rt) < 3, 
everywhere on C and, in particular, on K. Since (j) G C'^{K), we can write 

— log(l + 2rt)(9(90(;z) < / C-dd(p{z), 
K^^ Jk 

= C||-990(^)||ii(^) 

TT 

where C is independent of A^, K, and z, including z on the real line, and the norm 
\\^dd(t){z)\\Li(^K) depends only on K. So then we have that 

{En{z) dx A dy, (p^z)) < Ck 
where Ck is a constant which depends only on K. We now have want we want: 

^{En{z) dx A dy, (t){z)) < ^Ck = O(iV-i). 

Note that when we consider compact sets K that include part of the real line, the weak 
limit is the only result we have. We do not get a strong result because the derivatives of 
r, s, t, and therefore En{z) blow up near the real line. When we find the weak limit and move 
the dd from the log term to the term as we did above, we avoid this problem: only the 
derivatives of r and t blow up near the real line, not the values of the functions themselves. 



24 



BRIAN MACDONALD 



References 

[BBL92] E. Bogomolny, O. Bohigas, and P. Leboeuf. Distribution of roots of random polynomials. Physical 

Review Letters, 68(18):2726-2729, 1992. 

[BBL96] E. Bogomolny, O. Bohigas, and P. Leboeuf. Quantum chaotic dynamics and random polynomials. 
Journal of Statistical Physics, 85(5):639-679, 1996. 

[BSZOOa] P. Bleher, B. Shiffman, and S. Zclditch. Poincare-Lelong Approach to Universality and Scaling of 
Correlations Between Zeros. Communications in Mathematical Physics, 208(3):771-785, 2000. 

[BSZOOb] P. Bleher, B. Shiffman, and S. Zelditch. Universality and scaling of correlations between zeros on 
complex manifolds. Inventiones Mathematicae, 142(2):351-395, 2000. 

[DSZ04] M.R. Douglas, B. Shiffman, and S. Zelditch. Critical Points and Supersymmctric Vacua I. Com- 
munications in Mathematical Physics, 252(1) :325-358, 2004. 

[EK95] A. Edclman and E. Kostlan. How many zeros of a random polynomial are real? American Math- 
ematical Society, 32(1):1 37, 1995. 

[Han96] J.H. Hannay. Chaotic analytic zero points: exact statistics for those of a random spin state. J. 
Phys. A: Maih. Gen, 29(5):L101 L105, 1996. 

[IZ97] I. Ibragimov and O. Zcitoimi. On Roots of Random Polynomials. Transactions of the American 
Mathematical Society, 349(6):2427 2441, 1997. 

[Kac48] M. Kac. On the Average Number of Real Roots of a Random Algebraic Equation (II). Proceedings 
of the London Mathematical Society, 2(6):401, 1948. 

[Mac09] B. Macdonald. Density of Complex Critical Points of a Real Random Polynomial in Several Vari- 
ables, in final preparation, 2009. 

[Pro96] T. Prosen. Exact statistics of complex zeros for Gaussian random polynomials with real coefficients. 
Journal of Physics A: Mathematical and Ceneral, 29(15):4417-4423, 1996. 

[Ric54] S.O. Rice. Mathematical analysis of random noise. Selected Papers on Noise and Stochastic Pro- 
cesses, pages 133-294, 1954. 

[SV95] L.A. Shepp and R.J. Vanderbei. The Complex Zeros of Random Polynomials. Transactions of the 
American Mathematical Society, 347(ll):4365-4384, 1995. 

[SZ99] B. Shiffman and S. Zelditch. Distribution of Zeros of Random and Quantum Chaotic Sections of 
Positive Line Bundles. Communications in Mathematical Physics, 200(3) :661-683, 1999. 
E-mail address: binac@jhu.edu 



